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Abstract 

Let V be an 7i-polytope in R d with vertex set V. The vertex centroid is defined as the 
average of the vertices in V. We prove that computing the vertex centroid of an 7i-polytope 
is #P-hard. Moreover, we show that even just checking whether the vertex centroid lies 
in a given halfspace is already #P-hard for W-polytopes. We also consider the problem of 
approximating the vertex centroid by finding a point within an e distance from it and prove 
this problem to be #P-easy by showing that given an oracle for counting the number of ver- 
tices of an H-polytope, one can approximate the vertex centroid in polynomial time. We also 
show that any algorithm approximating the vertex centroid to any "sufficiently" non-trivial 
(for example constant) distance, can be used to construct a fully polynomial approximation 
scheme for approximating the centroid and also an output-sensitive polynomial algorithm 
for the Vertex Enumeration problem. Finally, we show that for unbounded polyhedra the 
vertex centroid can not be approximated to a distance of d^~ s for any fixed constant 5 > 0. 

1 Introduction 

An intersection of a finite number of closed halfspaces in M d defines a polyhedra. A polyhedra 
can also be represented as conv(V) + cone(Y), the Minkowski sum of the convex hull of a finite 
set of points V and the cone of a finite set of rays. A bounded polyhedra is called a polytope. 
In what follows, we will discuss mostly polytopes for simplicity and refer to the unbounded case 
explicitly only towards the end. 

Let V be an W-polytope in M with vertex set V. Various notions try to capture the essence 
of a "center" of a polytope. Perhaps the most popular notion is that of the center of gravity of 
V . Recently Rademacher proved that computing the center of gravity of a polytope is #P-hard 
[7]. The proof essentially relies on the fact that the center of gravity captures the volume of 
a polytope perfectly and that computing the volume of a polytope is #P-hard [3]. Note that, 
polynomial algorithms exist that approximate the volume of a polytope within any arbitrary 
factor pEj. It is also easy to see that the center of gravity can be approximated by simply 
sampling random points from the polytope, the number of samples depending polynomially on 
the desired approximation (See Algorithm 5.8 of [1]). 

In this paper we study a variant of the notion of "center" defined as the centroid (average) of 
the vertices of P. Despite being quite a natural feature of polytopes, this variant seems to have 
received very little attention both from theoretical and computational perspectives. Throughout 
this paper we will refer to the vertex centroid just as centroid. The reader should note that 
in popular literature the word centroid refers more commonly to the center of gravity. We 
nevertheless use the same terminology for simplicity of language. Our motivation for studying 
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the centroid stems from the fact that the centroid encodes the number of vertices of a polytope. 
As we will see, this also makes computing the centroid hard. 

The parallels between centroid and the center of gravity of a polytope mimic the parallels 
between the volume and the number of vertices of a polytope. Computing the volume and 
the number of vertices are both #P-complete ([21 13 [6]) and so are the problems of computing 
the corresponding centroids ([7], Theorem [T]). The volume can be approximated quite well but 
approximating the number of vertices of a polytope is an interesting open problem. Similarly, 
the center of gravity can be approximated quite well but (as we will see in this paper) obtaining 
a polynomial algorithm for approximating the centroid would be a very interesting achievement. 

The problem of enumerating vertices of an W-polytope has been studied for a long time. 
However, in spite of years of research it is neither known to be hard nor is there an output 
sensitive polynomial algorithm for it. A problem that is polynomially equivalent to the Vertex 
Enumeration problem is to decide if a given list of vertices of an W-polytope is complete pQ. In 
this paper we show that any algorithm that approximates the centroid of an arbitrary polytope 
to any "sufficiently" non-trivial distance can be used to obtain an output sensitive polynomial 
algorithm for the Vertex Enumeration problem. 

The main results of this paper are the following: 

• Computing the centroid of an H-polytope is #P-hard. 

• Even just deciding whether the centroid of an "H-polytope lies in a halfspace remains 
#P-hard. 

• Approximating the centroid of an "H-polytope is #P-easy. 

• Any algorithm approximating the centroid of an arbitrary polytope within a distance 
dz~ s can be used to obtain a fully polynomial approximation scheme for the centroid 
approximation problem and also an output sensitive polynomial algorithm for the Vertex 
Enumeration problem. 

• There is no polynomial algorithm that approximates the vertex centroid of arbitrary TL- 
polyhedron within a distance dz~~ s for any fixed constant 5 > 0, unless P = NP. 

We should remark that for the approximation of centroid, we only consider polytopes (and 
polyhedra) whose vertices lie inside a unit hypercube. To see how this assumption can easily 
be satisfied, notice that a halfspace h can be added to a polyhedron P such that P Pi h is 
bounded and the vertices of P are preserved in P fl h. Also, such a halfspace can be found 
in polynomial time from the inequalities defining P. Once we have a polytope in M. d , solving 
2d linear programs gives us the width along each coordinate axis. The polytope can be scaled 
by a factor depending on the width along each axis to obtain a polytope all whose vertices lie 
inside a unit hypercube. In case we started with a polyhedron P, the scaled counterpart of the 
halfspace h that was added can be thrown to get back a polyhedron that is a scaled version 
of P and all whose vertices lie inside the unit hypercube. In subsection 12.21 we provide further 
motivation for this assumption. 

Since all the vertices of the polytope (or polyhedron) lie inside a unit hypercube, picking 
any arbitrary point from inside this hypercube yield a d 2 -approximation of the vertex centroid. 
Thus, the last result should be contrasted to the fact that approximating the vertex centroid 
within a distance of is trivial. Also, even though we discuss only polytopes i.e. bounded 
polyhedra in subsections 12.11 and 12.21 the results and the proofs are valid for the unbounded 
case as well. We discuss the unbounded case explicitly only in subsection 12.31 
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2 Results 



2.1 Exact Computation of the Centroid 

The most natural computational question regarding the centroid of a polytope is whether we 
can compute the centroid efficiently. The problem is trivial if the input polytope is presented 
by its vertices. So we will assume that the polytope is presented by its facets. Perhaps not 
surprisingly, computing the centroid of an W-polytope turns out be #P-hard. We prove this by 
showing that computing the centroid of an "H-polytope amounts to counting the vertices of the 
same polytope, a problem known to be #P-hard. 

Theorem 1 Given an Ti-polytope V C R d , it is #P-hard to compute its centroid c(V). 

Proof. Embed V in K +1 by putting a copy of V in the hyperplane x^+i = 1 and making 
a pyramid with the base V and apex at the origin. Call this new polytope Q. Treating the 
direction of the positive x^+i-axis as up, it is easy to see that the centroid of the new polytope 
lies at a height 1 — ^j-j- iff the number of vertices of V is n. Thus any algorithm for computing 
the centroid can be run on Q and the number of vertices of V can be read off the (d + l)-st 
coordinate. □ 

Suppose, instead, that one does not want to compute the centroid exactly but is just in- 
terested in knowing whether the centroid lies to the left or to the right of a given arbitrary 
hyperplane. This problem turns out to be hard too, and it is not difficult to see why. 

Theorem 2 Given an Ti-polytope V C K d and a hyperplane h = {a ■ x = b}, it is #P-hard to 
decide whether a ■ c(V) < b. 

Proof. Consider the embedding and the direction pointing upwards as used in the proof of 
Theorem [TJ Given an oracle answering sidedness queries for the centroid and any arbitrary 
hyperplane, one can perform a binary search on the height of the centroid and locate the exact 
height. The number of queries needed is only logarithmic in the number of vertices of V . □ 



2.2 Approximation of the Centroid 

As stated before, even though computing the gravitational centroid of a polytope exactly is 
#P-hard, it can be approximated to any precision by random sampling. Now we consider 
the problem of similarly approximating the vertex centroid of an W-polytope. Let dist(x, y) 
denote the Euclidean distance between two points x,y £ M d . We are interested in the following 
problem: 

Input: "H-polytope Pct l! and a real number e > 0. 
Output: p el d such that dist(c(P),p) < e. 

We would like an algorithm for this problem that runs in time polynomial in the number of 
facets of P, the dimension d and ~. Clearly, such an algorithm would be very useful because 
if such an algorithm is found then it can be used to test whether a polytope described by m 
facets has more than n vertices, in time polynomial in m, n and the dimension d of the polytope 
by setting e < — ^y^j ■ This in turn would yield an algorithm that computes the number 
of vertices n of a d-dimensional polytope with m facets, in time polynomial in m, n and d. As 
stated before, a problem that is polynomially equivalent to the Vertex Enumeration problem is 
to decide if a given list of vertices of an "H-polytope is complete pQ. Clearly then, a polynomial 



3 



approximation scheme for the centroid problem would yield an output-sensitive polynomial 
algorithm for the Vertex Enumeration problem. 

Also, the problem of approximating the centroid is not so interesting if we allow polytopes 
that contain an arbitrarily large ball, since this would allow one to use an algorithm for ap- 
proximating the centroid with any guarantee to obtain another algorithm with an arbitrary 
guarantee by simply scaling the input polytope appropriately, running the given algorithm and 
scaling back. So we will assume that the polytope is contained in a unit hypercube in W 1 . 

Now we prove that the problem of approximating the centroid is #P-easy. We do this by 
showing that given an algorithm that computes the number of vertices of an arbitrary polytope 
(a #P-complete problem), one can compute the centroid to any desired precision by making 
a polynomial (in -, the number of facets and the dimension of the polytope) number of calls 
to this oracle. Notice that in the approximation problem at hand, we are required to find a 
point within a d-ball centered at the centroid of the polytope and radius e. We first modify 
the problem a bit by requiring to report a point that lies inside a hypercube, of side length 
2e, centered at the centroid of the polytope. (The hypercube has a clearly defined center of 
symmetry, namely its own vertex centroid.) To see why this does not essentially change the 
problem, note that the unit hypercube fits completely inside a (i-ball with the same center and 
radius We will call any point that is a valid output to this approximation problem, an 
e-approximation of the centroid c(P). 

Given an W-polytope P and a hyperplane {a-x = b} that intersects P in the relative interior 
and does not contain any vertex of P, define Pi and P2 as follows: 

P 1 = Pn{x\a-x<b} 
P 2 = Pn{x\a-x>b} 

Let V\ be the common vertices of Pi and P, and V2 be common vertices of P2 and P. 
The following lemma gives a way to obtain the e-approximation of the centroid of P from the 
e-approximations of the centroids of V\ and V2 ■ 

Lemma 1 Given P, V±, V2 defined as above, let n\ and 712 be the number of vertices in V\ and 
V2 respectively. If c\ and c 2 are e-approximations of the centroids of V\ and V2 respectively, 
then c = Wl ^^^ C2 is an e-approximation of the centroid c* of P. 

Proof. Let Cy be the j-th coordinate of q for i £ {1, 2}. Also, let c* be the actual centroid of 
Vi with c*j denoting the j'-th coordinate of c*. Since Cj approximates c* within a hypercube of 
side- length 2e, for each j € {1, ■ ■ ■ , d} we have 



Cjj e < Cij < c,ij + e 

Also, since c* is the centroid of P, 

„ _ nic* + n 2 c* 2 



ni + n 2 



Hence, for each coordinate c* of c* we have 



n 1 (c 1] -e)+n 2 (c 2j ^e) < ^ < n x (cij + e) + n 2 (c 2j + e) 

ni+n 2 - 3 - ni+n2 

n\c\j+n 2 c 2 j _ < c * < re l C lj + n 2 c 2j 

n 1+ n 2 - - j - ni + ri2 

Cj -e < c* < Cj + e 

c* - e < Cj < c* + e 
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□ 



Now to obtain an approximation of the centroid, we first slice the input polytope P from 
left to right into \ slices each of thickness at most e. Using standard perturbation techniques 
we can ensure that any vertex of the input polytope does not lie on the left or right boundary 
of any slice. For each slice any point in the interior gives us an e-approximation of the vertices 
of P that are contained in that slice. We can compute the number of vertices of P lying in this 
slice by using the oracle for vertex computation and then using the previous Lemma we can 
obtain the centroid of P. Thus we have the following theorem: 

Theorem 3 Given a polytope P contained in the unit hypercube, the e-approximation of the 
centroid of P can be computed by making a polynomial number of calls to an oracle for computing 
the number of vertices of a polytope. 

Now we present a bootstrapping theorem indicating that any "sufficiently" non-trivial ap- 
proximation of the centroid can be used to obtain arbitrary approximations. For the notion 
of approximation let us revert back to the Euclidean distance function. Thus, any point x 
approximating the centroid c within a parameter e satisfies dist(x, c) < e. As before we assume 
that the polytope V is contained in the unit hypercube. Since the polytope is thus contained 
in a hyperball with origin as its center and radius at most , any point inside V approximates 
the centroid within a factor \fd. Before we make precise our notion of "sufficiently" non-trivial 
and present the bootstrapping theorem, some preliminaries are in order. 

Lemma 2 Suppose (x,y), (u,u) G R 2d , where x,y,u G R d , then 

ii X + V n <r - ( x >y)\\ 

ll« 2 II < ^ 
where \ \ ■ \ \ is the Euclidean norm. 

The proof of the above lemma is easy and elementary, and hence we omit it here. Next, 
consider the product of two polytopes. Given ci-dimensional polytopes V, Q the product V x Q 
is defined as the set {(x,y)\x 6 P,j/ 6 Q}. The facet defining inequalities of the product of 
P, Q can be computed easily from the inequalities defining P and Q. 

P = {x|^ix<6i}, 

Q = {y\A 2 y < b 2 } 
=> PxQ = {(x, y)\A lX < h,A 2 y < b 2 }, 



where A 1 G R miXdl ,A 2 G R m ^,x G R dl ,y G M. d2 ,h G R miXl ,b 2 G M m2Xl . 

It is easy to see that the number of vertices of V x Q is the product of the number of vertices 
of V and that of Q, and the number of facets of V x Q is the sum of the number of facets of V 
and that of Q. Moreover, the dimension of V x Q is the sum of the dimensions of V and that 
of Q. 

Observation 1 If c is the centroid of a polytope P then (c, c) is the centroid of P x P. 

Suppose we are given an algorithm for finding e-approximation of an arbitrary polytope 
contained in the unit hypercube. For example, for the simple algorithm that returns an arbitrary 
point inside the polytope, the approximation guarantee is y ^ L . We consider similar algorithms 
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whose approximation guarantee is a function of the ambient dimension of the polytope. Now 
suppose that for the given algorithm the approximation guarantee is f{d). For some parameter 

k times 

' * 1 k 

k consider the /c-fold product of P with itself P x ■ ■ • x P, denoted by P . Using the given 
algorithm one can find the f(2 k d) approximation of P k and using Lemma [2] one can then find 
the y k -approximation of P. This gives us the following bootstrapping theorem: 

Theorem 4 Suppose we are given an algorithm that computes a approximation for any 
polytope contained in the unit hypercube in polynomial time, where g(.) is an unbounded mono- 
tonically increasing function. Then, one can compute an e- approximation in time polynomial 
in the size of the polytope and g~ 1 (^r). 

In particular, if we have an algorithm with d^~ s approximation guarantee for finding the 
centroid of any polytope for some fixed constant 5 > 0, then this algorithm can be used to 
construct a fully polynomial approximation scheme for the general problem. 

2.3 Approximating centroid of a polyhedron is hard 

The reader should note that the analysis of subsections 12.11 and 12.21 remains valid even for 
the unbounded case (polyhedra). Even though we do not have any idea about the complexity 
of approximating the centroid of a polytope, now we show that for an arbitrary unbounded 
polyhedron the vertex centroid can not be ^"^-approximated for any fixed constant 5 > 0. To 
show this we first prove that for an 7Y-polyhedron P C M. d the vertex centroid of P can not be 
^-approximated in polynomial time unless P = NP. This together with Theorem 0] completes 

the proof for hardness of d 2 ~ s - approximation of the centroid of an W-polyhedron. 

The proof sketch is as follows: Given a boolean CNF formula <j>, we construct a graph G((j>) 
such that G((j>) has a "long" negative cycle if and only if <j) is satisfiable. For a given graph 
G we define a polyhedron P(G) such that every negative cycle in G is a vertex of P(G) and 
vice- versa. From the properties of the vertex centroid of this class of polyhedron, we then prove 
that for any formula </>, ^-approximating the vertex centroid of P(G((f>)) would reveal whether 
4> is satisfiable or not. 

Before we proceed further, we would like to remark that for the graph corresponding to a 
given 3-CNF formula, we use the same construction that appears in [5] except for different edge 
weights. We nevertheless include the discussion here for completeness. Our main contribution 
is the negative-flow polyhedron associated with a directed graph and a set of edge weights and 
the analysis of the vertices of this polyhedron leading to Theorem [5j 

2.3.1 Graph of a CNF formula 

Given a 3-CNF 4>{x\, . . . ,xn) = C\ A . . . A C m , we construct a weighted directed graph G = 
(V,E), on |V| = 5£Jli \Cj\+m-n+l = 16m-n+l vertices and \E\ = 6J2T=i \ C j\ +1 = 1 & m + 1 
arcs (where \Cj\ denotes the number of literals appearing in clause Cj, i.e. 3 for a 3-CNF), as 
follows. For each literal £ = P appearing in clause Cj, we introduce two paths of three arcs 
each: V(£) = {p(£) , a(£) , b(£) , q{£)) , and V (£) = (r(£),b' (£), a' (£), s(£)). The weights of these 
arcs are set as follows: 

w((p(£), a(£))) = ±, w((a(£), b(£))) = -\, w{(b(£),q(£))) = 0, 
w((r(£),b'(£)))=0, w {(b\£),a'{£))) = -\, w((b> (£), s(£))) = \. 
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Figure 1: An example of the graph construction for CNF C = {x\ V x 2 V 23) A [x\ V X2 V X3) A 
(xi V x 2 Vx 3 ). 

These edges are connected in G as follows (see Figure [T] for an example) : 

G = v g 1 vi Q 2 v 2 . ■ ■ v n -i Q n v n g[ v[ g' 2 v' 2 ... v' m _ 1 Q' m v' m , 

where vo,vi, . . . , v n , v' m _ 1 ,v' m are distinct vertices, each ft = J'i V Zi, for i = 1, . . . , n, 

consists of two parallel chains 3^ = AjV(xl) and i% = AjV(xl) between V{-\ and Uj, and each 
IC ■ I 

0'. = v^T'(^), for j = 1, . . . , to, where . . . are the literals appearing in Cj. 

Finally we add the arc (v' m ,Vo) with weight —1, and identify the pairs of nodes {a(£),a'(£)} 
and {b(£),b'(£)} for all £, (i.e. a(£) = a{£') and 6(^) = 6(f) define the same nodes). 

We will argue now that all negative cycles of G have weight —1. Clearly the arcs (a(£), b(£)) 
and (b'(£), a! {£)) form a directed cycle of total weight —1, for every literal occurrence £. There 
are YljLi = ^ m such cycles, corresponding to a subset X C V(P(G, w)). 

Call a cycle of G long if it contains the vertices vo,Vi,...,v n ,v^,..., v' m __ 1 ,v' m . Any long 
cycle has weight —1. The crucial observation is the following. 

Claim 1 Any negative cycle C £ C~(G,w) \ X must be long. 

Proof. Consider any cycle C X, and let us write the traces of the nodes visited on the 
cycle (dropping the literals, and considering a, a' and b,b' as different copies), without loss of 
generality as follows: 

pabpabp---aasb'asb •••bbpab---p. 

Note that the sequences a' a and b b' are not allowed since otherwise C contains a cycle from 
X. 

Let us compute the distance (i.e., the total weight) of each node on this sequence starting 
from the initial p. Call the subsequences 0! a' and b b' , a- and 6-jumps respectively. Then it 
is easy to verify that each a-jump causes the distance to eventually increase by 1 while each 
6-jump keeps the distance at its value. More precisely, the distance at a node x in the sequence 
is given by d(x) = t + do — 5(x), where t is the number of a-jumps appearing upto x, and 

'0 ii x £ {p, s}, 

. _ 5 if x = a, g, , _ f 1 if arc (v' m , vo) appears on the path from p to x 

— 5 ifx = o', \0 otherwise. 

if £ = & = &', 
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One also observes that, if the sequence has a 6-jump, then it must also contain an o-jump. Thus 
it follows from the definition of d{x) that any cycle with a jump must be non-negative. So the 
only possible negative cycle not in X must be long. □ 

If we are interested only in the negative cycles then it is clear that long cycles are exactly 
those negative cycles that contain the edge (v' m ,vo). In what follows this edge will be called 
the marked edge, the negative cycles containing this marked edge will be called long negative 
cycles and the other negative cycles will be called short negative cycles. 

By Claim [TJ checking of V(P(G, w)) = X is equivalent to checking if G has a long cycle. It 
is easy to see that the latter condition is equivalent to the non-satisfiability of the input CNF 
formula (j) ( see [S] f° r more details). 

To summarize this subsection, for a 3-CNF with m clauses we can obtain a graph G with 
18m + 1 edges such that G has 3m short negative cycles and finding whether G has any negative 
cycle (called long negative cycle) containing one special edge is NP-complete. Furthermore, each 
of these negative cycles has total weight —1. 



2.3.2 The polyhedron of negative-weight flows of a graph 



Given a directed graph G = (V, E) and a weight function w : E 
following polyhedron: 



P(G,w) 



y e 



Uvu 

v:(u,v)€E v:(v,u)gE 



on its arcs, consider the 



V ueV 



(N) 



Vuv > 



V (u,v) G E 



If we think of w UtV as the cost/profit paid for edge (u, v) per unit of flow, then each point of 
P(G,w) represents a negative-weight circulation in G, i.e., assigns a non-negative flow on the 
arcs, obeying the conservation of flow at each node of G, and such that total weight of the flow 
is strictly negative. 

For a subset X C E, and a weight function w : E i— ► R, we denote by w(X) = X^ex^e, 
the total weight of X. For X C E, we denote by xPO £ {Oi the characteristic vector of 
X: Xe(X) = 1 if and only if e G X, for e G E. The following theorem states that the vertex set 
V(P(G, w)) of P(G, w) is in one-to-one correspondence with the negative cycles of the graph G. 



Theorem 5 

Then 



Let G = (V, E) be a directed graph and w : E 



be a real weight on the arcs. 



V(P(G,w)) 



w(C) 



X (C): CeC~(G,w) 



(1) 



We prove this theorem by proving a series of claims about the graph G. Let m = \E\ and 
n = \V\. It is easy to verify that any element y G M s of the set on the right-hand side of 
([1]) belongs to P(G,w). Moreover, any such x = —x(C)/w(C), for a cycle C, is a vertex of 
P(G, w) since there are m linearly independent inequalities of P(G, w) tight at x, namely: the 
conservation of flow equations (F) at \C\ — 1 vertices of C, the equation X^eeC w eVe = — 1> an< ^ 
m — \C\ equations y e = 0, for e G E \ C. 

To prove the opposite direction, let y G M. E be a vertex of P(G, w). Let Y = {e G E : y e > 
0}. The proof follows from the following 3 claims. 
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Claim 2 The graph (V, Y) is the disjoint union of strongly connected components. 



Proof. Consider an arbitrary strongly connected component X in this graph, and let X~ 
be the set of components reachable from X (including X). Summing the conservation of flow 
equations corresponding to all the nodes in X~ implies that all arcs going out of X~ have a 
flow of zero. □ 



Claim 3 There exists no cycle C G C°(G,w) such that CCF. 
Proof. If such a C exists, we define two points y' and y" as follows. 



/ 



y e + e, if e G C „ ( y e - e, if e G C 



y = < y 

y e , otherwise, ' ' [ y e , otherwise, 

for some sufficiently small e > 0. Then y',y" clearly satisfy (F). Moreover, (N) is satisfied by y' 
since 

^2 Wey ' e = ^2 WeVe + ^2 We ( ye + e ) = WeVe + w ( c ) e = ~ L 

e£E e^C eeC e£E 

Similarly for y" . Thus y',y" G P(G,w) and y = (y' + y")/2 contradicting that y is a vertex. □ 

Claim 4 There exist no distinct cycles C\,Ci G C~(G, w) L)C + (G,w) such that C\ U Ci Q Y. 
Proof. If such C\ and C2 exist, we define two points y' and y" as follows. 



y' e 



y e + ei, ifeGCi\C 2 

y e + e 2 , ifeGC 2 \Ci , 

y e + ei + e 2 , if e G C\ n C 2 2/6 

y e , otherwise, 

w(Ci) ( 



y e -ei, ifeGCi\C 2 

y e - e 2 , if e G C 2 \ C\ 

y e -ei-e 2 , if e G C\ n C 2 

y e , otherwise, 



where ei = — ^|^| e 2 , for some sufficiently small e 2 > (in particular, to insure non-negativity 

of y' ,y" , e 2 must be upper bounded by the minimum of min{y e : e G C 2 \ Ci}, |^|^|| min{y e : 

e G Ci \ C 2 }, and |^^^Lyy min{y e : e G C\ n C 2 }). Then it is easy to verify that y',y" 
satisfy (F). Moreover, (N) is satisfied with y'since 

^2w e y' e = ^2 w eVe+ ^2 w e{Ve + ei)+ ^ w e {y e + e 2 ) 

e&E e£C 1 UC 2 e£d\C 2 eeC 2 \Ci 

+ ^2 w e(Ue + ^1 + ^2) =/ J WeVe + w{C\)ex + w(C 2 )e 2 = -1. 
eednC 2 eGE 

Similarly for y". Thus 2/,y" G P(G,w) and y = (y' + y")/2 contradicting that y is a vertex of 
P(G,w). □ 

The above 3 claims imply that the graph (V, 1") consists of a single cycle C and a set of 
isolated vertices V \ V(C). Thus y e = for e C. By (F) we get that y e is the same for all 
e G C, and by (N) we get that y e = —l/w(C) for all e G C, and in particular that C G C~ (G, w). 
This completes the proof of Theorem [5j 

Since, in the graph G arising from a 3-CNF formula, every negative cycle has weight exactly 
— 1, Theorem [5] implies that the vertices of P(G,w) are exactly the characteristic vectors of the 
negative cycles of G. Recall that finding whether G has any long negative cycle, i.e. a negative 
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cycle containing a marked edge (say e), is NP-complete. Further recall that for a 3-CNF formula 
with m clauses the constructed graph G has 18m + 1 edges and 3m trivial short negative cycles. 
Consequently, the polyhedron P(G, w) that is finally obtained has dimension 18m + 1 and 3m 
trivial vertices corresponding to the short negative cycles of G. 

Now, if there are no long negative cycles then the vertex centroid of P(G, w) has value in 
the coordinate corresponding to the edge e. For simplicity, we will refer to this coordinate axis 
as x e . On the other hand, if there are K > 1 long negative cycles in G then in the centroid 
Xe = K+3m — 3m+i ' This implies that having an e- approximation for the centroid of P(G,w) 
f° r 6 ^ 2(3m+i) wou ld reveal whether or not P(G, w) has a non-trivial vertex and hence whether 
or not G has a long negative cycle. Thus we have the following theorem: 

Theorem 6 There is no polynomial algorithm that computes a ^-approximation of the vertex 
centroid of an arbitrary 7i -polyhedron P C M d , unless P = NP. 

An immediate consequence of Theorem H] and Theorem [6] is that there is no polynomial 
algorithm that computes any "sufficiently non-trivial" approximation of the vertex centroid of 
an arbitrary W-polyhedron unless P = NP. More formally, 

Corollary 1 There is no polynomial algorithm that d?~ s -approximates the centroid of an ar- 
bitrary d- dimensional TC -polyhedron for any fixed constant 5 > unless P = NP. 

3 Open Problems 

Although we can show that for unbounded polyhedra almost any non-trivial approximation of 
the vertex centroid is hard, we can not make a similar statement for the bounded case (i.e. 
polytopes. One interesting variant of Theorem [2] would be to consider a ball of radius r instead 
of a halfspace. If containment of vertex centroid in a ball of radius r can be decided in time 
polynomial in the number of inequalities defining the polytope, the dimension and r then one can 
perform a sort of random walk inside the polytope and approximate the centroid in polynomial 
time. We leave out the details of this random walk since we do not have a method to check 
containment inside a ball. 

References 

[1] D. Avis, D. Bremner, and R. Seidel. How good are convex hull algorithms? Corn-put. Geom., 
7:265-301, 1997. 

[2] M. E. Dyer. The complexity of vertex enumeration methods. Mathematics of Operations Research, 
8(3):381-402, 1983. 

[3] M. E. Dyer and A. M. Frieze. On the complexity of computing the volume of a polyhedron. SIAM 
J. Comput., 17(5):967-974, 1988. 

[4] R. Kannan, L. Lovasz, and M. Simonovits. Random walks and an o*(n 5 ) volume algorithm for convex 
bodies. Random Structures and Algorithms, 11(1):1— 50, December 1998. 

[5] L. Khachiyan, E. Boros, K. Borys, K. M. Elbassioni, and V. Gurvich. Generating all vertices of a 
polyhedron is hard. In SODA, pages 758-765. ACM Press, 2006. 

[6] N. Linial. Hard enumeration problems in geometry and combinatorics. SIAM J. Algebraic Discrete 
Methods, 7(2):331-335, 1986. 

[7] L. Rademacher. Approximating the centroid is hard. In Symposium on Computational Geometry, 
pages 302-305, 2007. 



10 



